Mapping spacetimes with LISA: inspiral of a test-body in a 'quasi-Kerr' field. 
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The future LISA detector will constitute the prime instrument for high-precision gravitational 
wave observations. Among other goals, LISA is expected to materialise a 'spacetime-mapping' 
program, that is to provide information for the properties of spacetime in the vicinity of supermassive 
black holes which reside in the majority of galactic nuclei. Such black holes can capture stellar-mass 
compact objects, which afterwards slowly inspiral under the emission of gravitational radiation. 
The small body's orbital motion and the associated waveform observed at infinity carry information 
about the spacetime metric of the massive black hole, and in principle it is possible to extract 
this information and experimentally identify (or not!) a Kerr black hole. In this paper we lay the 
I/-) ■ foundations for a practical spacetime-mapping framework. Our work is based on the assumption 

that the massive body is not necessarily a Kerr black hole, and that the vacuum exterior spacetime is 
stationary-axisymmetric, described by a metric which deviates slightly from the known Kerr metric. 
£NJ ■ We first provide a simple recipe for building such a 'quasi-Kerr' metric by adding to the Kerr metric 

the leading order deviation which appears in the value of the spacetime's quadrupole moment. We 
then study geodesic motion of a test-body in this metric, mainly focusing on equatorial orbits, but 
■ also providing equations describing generic orbits formulated by means of canonical perturbation 

theory techniques. We proceed by computing approximate 'kludge' gravitational waveforms which 
we compare with their Kerr counterparts. We find that a modest deviation from the Kerr metric is 
sufficient for producing a significant mismatch between the waveforms, provided we fix the orbital 
parameters. This result suggests that an attempt to use Kerr waveform templates for studying 
extreme mass ratio inspirals around a non-Kerr object might result in serious loss of signal-to-noise 
ratio and total number of detected events. The waveform comparisons also unveil a 'confusion' 
problem, that is the possibility of matching a true non-Kerr waveform with a Kerr template of 
' different orbital parameters. 
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V/-) ' I. INTRODUCTION 

o ' 

During the past several years astronomical observations have provided almost indisputable evidence in favour of 
U the existence of 'dark' supermassive objects (with a mass spectrum ~ 10 s — 10 9 M Q ) in the majority of galactic cores, 
including our own Milky Way Q. Conventional wisdom dictates that these objects should be Kerr black holes as 
&)[)■ described by General Relativity. This belief is (so far) rather based on our faith in General Relativity itself rather 
than on hard 'experimental' evidence 1 . Other candidate massive objects have been proposed (such as soliton stars, 
boson stars, gravastars and P-stars but these are treated with scepticism by the community as they involve 'exotic' 
physics. 

The future LISA gravitational wave detector Q is expected to be able to give a definitive answer to whether these 
massive objects are Kerr black holes or not. It is commonly accepted [f| that LISA will have the potential of 'mapping' 
the spacetime in these objects' close vicinity. This can be achieved by detection of gravitational radiation emitted 
during the inspiral of a stellar-mass compact body around the supermassive object. It is expected that LISA will be 
capable of detecting from several to about a thousand 6] of these extreme mass ratio inspirals (EMRIs), during a 
3 — 5 years mission. For the actual detection and subsequent extraction of the system's physical properties, LISA's 
data analysis will rely heavily on matched filtering. This method is based on cross-correlating the detector's noisy 
output with a pre-selected bank of waveform templates which should be accurate representations of the true signal. 
Among anticipated gravitational wave sources, EMRIs are thought of as one of the 'cleanest': they can be accurately 
modelled as a binary system of a 'test-body' orbiting a Kerr black hole, evolving under its own gravitational dynamics, 
without receiving any other significant 'environmental' influence 2 . 
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1 Note however that the existence of an event horizon seems to be required in order to explain the spectral properties of X-ray emission 
by Galactic black hole binaries and AGNs, see for example Ref. 0] 

2 This should be true for the majority of galactic supermassive black holes which are in a "quiescent" accretion state. For the small 
fraction of strongly accreting black holes (like the ones in AGNs) the presence of an accretion disk may have non-negligible effects on a 
orbiting body (see for example Q) 
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The basic tool for studying such a system is black hole perturbation theory: the small body of mass (J, ~ 1 — 10M Q 
perturbing the spacetime of the much more massive black hole (of mass M ~ 10 6 ). For recent detailed reviews on 
EMRIs and further references on the subject we refer the reader to Refs. Q, 

In their great majority, EMRI-related studies take for granted the Kerr identity of the supermassive object. However, 
if we ever hope to materialise any spacetime-mapping program with LISA we should be prepared to dismiss (at least 
to some degree) this assumption. In fact, any claim that LISA will be able to probe the Kerr metric would require 
the use of realistic waveform templates which contain -in some suitable parametric form- deviations from the Kerr 
metric. The task of a LISA data analyst would be to quantify these deviations, most likely pointing to a null result. 

This idea was first put forward by Fintan Ryan in the late 90s. In a series of papers Ryan demonstrated 
how LISA could in principle construct a map of the massive body's spacetime by means of EMRIs observations. His 
method is based on writing the general stationary axisymmetric vacuum metric in terms of multipole moments as 
prescribed in Ref . ^l| • These can be mass moments Mg and current moments Si , both labelled by the angular integer 
eigenvalue i > 0. Given the moments, the metric takes the symbolic form [To|. Ill- 
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9ab ~ ^ r 2n+l {^2n^W + (M(<2ni S e<2n )} , ^ 2 „ +3 { S 2n+1 V% (6>) + (M e<2n +1 , S e<2n+1 ) } 



n=0 L 



(1.1) 



where V^ 2 (()) are known angular functions and angular brackets stand for the contribution of lower moments (if they 
exist for given n). 

These multipole moments will be clearly encoded in the geodesic equations of motion 3 for a test-body and conse- 
quently in the phase and amplitude of the emitted waveform. LISA (as any other gravitational wave detector) will 
be capable of tracking this phase with high precession and extract valuable information concerning the spacetime's 
multipolar structure. 

The key feature that makes this method so attractive is that Kerr spacetime is special, in the sense that all of its 
higher moments are 'locked' to the first two, the mass M = M and spin J = Si [U\. 

M t + iS t = M{iaf (1.2) 

where a = J/M is the familiar Kerr spin parameter. This relation is unique to Kerr holes and is a mathematical 
expression of the famous 'no-hair' theorem (l2l |. In practise this means that extracting just the three lowest moments, 
M, J and (say) the quadrupole moment M 2 , is sufficient for verifying (or disproving!) that the central body is indeed 
a Kerr black hole. Measurement of additional higher moments could pinpoint the identity of the body 4 . We should 
emphasise that a possible non-Kerr multipolar structure does not necessarily imply a non-Kerr identity for the massive 
object. A true Kerr hole with a substantial amount of material in its vicinity (say, in the form of an accretion disk) 
could effectively behave as a non-Kerr object. This scenario is not included in our present discussion but certainly 
merits a future detailed investigation. 

Ryan's multipole expansion scheme is certainly an elegant and powerful method for 'mapping spacetimes', never- 
theless there are some serious issues regarding its practicality. Computing gravitational waveforms for EMRIs requires 
(i) precise knowledge of a test-body's geodesic motion and (ii) a wave-emission formalism that allows the calculation 
of waveforms and fluxes, once the orbital motion is prescribed. For the case of a Kerr black hole, geodesic motion is 
well known and extensively studied |15| , while wave dynamics is studied with the help of the celebrated Teukolsky 
formalism p"ri| which encapsulates the dynamics of gravitational perturbations in one single master equation. How- 
ever, it is not always appreciated that the study of both orbital and wave dynamics is greatly simplified due to the 
speciality of Kerr spacetime. In particular, one is able to first decouple and then separate the dynamical equations 
thanks to the Petrov type-D character of the Kerr spacetime and a suitable choice of coordinate frame 

Life is considerably more complicated if we consider a non-Kerr spacetime with arbitrary multipolar structure. To 
begin with, according to (|l.lf> an accurate representation of the metric in the strong field regime (the most relevant 
part of the inspiral for LISA) requires the inclusion of a large number of multipoles. This point can be made clearer if 
we notice that the metric (|1.1|) is an expansion in 1/r around the flat Minkowski spacetime. A more suitable expansion 
would be around (say) the Schwarzschild spacetime. We will return to this point in future work. 

Complications first appear at the level of geodesic motion. The spacetime's stationary-axisymmetric character 
guarantees the existence of two integrals of motion, the energy E and angular momentum (along the symmetry axis) 



3 Multipolar expansion of the gravitational potential is a powerful technique in Newtonian gravity for performing geodesic measurements: 
the motion of a satellite around the Earth p rov ides information on our planet's 'bumpiness'. 

4 For example, a rotating massive Boson star has three independent moments and its higher moments certainly do not obey relation 
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L z . However, the third integral of motion (the famous Carter constant for the Kerr metric) will typically be lost, as 
the Hamilton- Jacobi equation is not fully separable (see ^3] f° r the solution of this equation in the Kerr case). Two 
integrals of motion are sufficient (together with initial conditions) for fully describing motion only for the special case 
of equatorial orbits. For generic orbits (which is the case relevant for LISA) the insufficient number of integrals of 
motion would force us to use the second-order geodesic equation of motion. This would result in a more complicated 
description of geodesic motion (as compared to the Kerr case), but technical complication is not the sole problem 
here. Fewer integrals of motion could easily translate into chaotic dynamical behaviour which in turn would have a 
great impact on LISA data analysis. 

The picture appears even more blurred when one attempts to study wave dynamics in a non-Kerr metric. Lack 
of decoupling for the perturbation equations inhibits the formulation of a 'Teukolsky-like' equation - the most ver- 
satile method in black hole perturbation theory is now lost. An alternative way of computing rigorous gravitational 
waveforms and fluxes would be the direct numerical time evolution of metric perturbations (after separation of the 
azimuthal <j) coordinate) which come in a package of ten coupled partial differential equations (combined with four 
gauge-condition equations) (l2^ . 

Motivated (and partially discouraged!) by the above complications, we adopt a somewhat different (albeit less 
general) approach in the present work. Based on the belief that the massive objects in galactic nuclei are most likely 
Kerr black holes, we only address the question: is the spacetime around these objects described by the exact Kerr 
metric, or by a slightly different metric ? We do not attempt to identify the source of this metric, in fact we are 
agnostic on this issue. Our less ambitious goal is to quantify the deviation from the Kerr metric. 

This point should be emphasised because none of the other candidate objects can be considered in any sense as 
being 'almost' a Kerr black hole. The idea of building a 'quasi-black hole' spacetime has been recently advocated by 
Collins & Hughes [l^. These authors constructed a 'bumpy' Schwarzschild black hole by adding a certain amount 
of quadrupole moment (in the form of a given mass distribution outside the black hole) and then studied equatorial 
orbits in the resulting spacetime. 

In our scheme we modify the quadrupole moment of the exact Kerr metric (this is the lowest moment where Kerr 
can be distinguished from the spacetime of another axisymmetric/stationary body). This procedure is performed 
in a natural way with the help of the well-known Hartle-Thorne exterior metric |2l|. The resulting 'quasi-Kerr' 
metric is customised for use in strong gravity situations and has the nice feature of smoothly reducing to the familiar 
Kerr metric in Boyer-Lindquist (BL) coordinates. Being arbitrarily close to the BL Kerr metric (which allows full 
separability of the Hamilton-Jacoby and wave equations) allows us to formulate first-order equations of motion and 
integrals of motion by perturbing the well known Kerr equations. As yet, we have not explored whether the same 
trick would lead to a 'quasi- Teukolsky' perturbation equation. For the purposes of the present investigation it is quite 
sufficient to use the so-called hybrid approximation H^jH^I and compute 'kludge' waveforms. 

The paper is structured as follows. In Section^we provide the prescription for building our quasi-Kerr metric. This 
is the paper's main result. Section IlIII is devoted to the solution of the Hamilton- Jacobi equation in the quasi-Kerr 
metric. We first discuss the issue of separation of variables for this equation and then proceed with the derivation 
of equations of motion based on standard canonical perturbation theory. Section IIII CI is dedicated to the study of 
equatorial orbits and comparison of orbital frequencies and periastron advance for Kerr and quasi-Kerr orbits. In 
Section llVI we compute approximate waveforms from test-bodies in equatorial orbits of the quasi-Kerr metric. These 
are compared to their Kerr analogues by calculating their overlap function in the LISA sensitivity band. In the 
same Section we touch upon the problem of 'confusion' between Kerr and quasi-Kerr waveforms. A summary and 
concluding discussion can be found in Section[V] Appendices containing some technical details, including a review of 
the action-angle formalism in the Kerr spacetime can be found at the end of the paper. Throughout the paper we 
use geometrised units G = c = 1 and adopt a { — , +, +, +} signature for the metric. We repeatedly use the labels 'K' 
and 'qK' for 'Kerr' and 'quasi-Kerr', respectively. 



As we mentioned in the Introduction, the central idea of our work is to replace the general multipolar expansion 
of an axisymmetric-stationary vacuum spacetime with a less general 'quasi-Kerr' metric. In the language of 
multipolc moments this means, 



II. BUILDING A QUASI-KERR SPACETIME 



Mi = Mf + 5M e , £>2 
S e = Sf + SSe, £>3 



(2.1) 
(2.2) 



where SM^, 8Si are assumed small. With respect to the first two leading moments M, J the spacetimes outside a Kerr 
hole and any other stationary axisymmctric rotating body are indistinguishable. Essentially, M and J can take the 



4 



same value for both bodies (the spin is presumably limited, for a Kerr hole a < M). This degeneracy is broken as 
soon as we take into account the next most important moment, the mass quadrupole M 2 . In our scheme we choose 
to take into account only the deviation in the quadrupole moment Q = M 2 and neglect any deviation in all higher 
moments. We introduce the dimensionlcss deviation parameter e as 5 



Q=Q K ~eM 3 , where Q K = (2.3) 

ivi 

Our quasi-Kerr metric takes the form, 

9»p = 9% + th af} + 0(SM £ > 4 , SS e > 3 ) (2.4) 

where g^a is the exact Kerr metric. Our goal is to determine the metric functions h a p(r, 9). 

This can be achieved by using the Hartle-Thorne (H-T) metric which describes the spacetime outside any slowly 
rotating axisymmetric and stationary body. Indeed, this metric is fully accurate up to the quadrupole moment (which 
scales quadratically with the body's spin), and it includes as a special case the Kerr metric (at 0(a 2 ) order). It is 
then possible to isolate the leading-order quadrupole moment deviation and deduce the h a p piece in (|2.4|l . Implicit in 
this prescription is the assumption that the massive body's exterior spacetime has nonzero higher moments induced 
by the body's rotation, in the same sense that a black hole has nonzero moments (with I > 0) only when a ^ 0. As 
a consequence, terms of order 0(ea,e 2 ) are beyond the desired accuracy and therefore neglected. Unlike the general 
multipole expansion our prescription by construction offers no information (regarding possible deviations from 

the Kerr values) on moments higher than the quadrupole. On the other hand our metric is not an expansion in inverse 
powers of r , instead is written in a compact fully relativistic form suitable for strong-field conditions. 

We move on to the calculation of h a b. For our purposes we shall only need the exterior H-T metric |2lj | which is 
expressed in terms of the dimensionless parameters, 

Then at 0{j 2 ) accuracy, the metric is given by (these expressions are taken from Ref. [24| after correcting a typo- 
graphical sign error in the F\ function), 

9u T = -(l-2M/r)[l+j 2 F 1+ qF 2 ] 

5j H r T = (l-2M/ry 1 [l+j 2 G 1 -qF 2 ] 

gfj = r 2 [l +] 2 H l+ qH 2 ] 

g*J = r 2 sin 2 9 [1+ f H 1 +qH 2 ] 

4 = ^sin 2 (2.6) 

The functions Fi_ 2 (r,9), Hi >2 (r,9) appearing in these expressions are given in Appendix 1X1 Note that the metric 
(12.01) is written in the coordinate frame {t, r, 9, (ft} originally used by Hartle & Thorne. Without loss of generality, we 
can write 

q - f + e (2.7) 
Clearly, the limit e — > corresponds to the H-T representation of Kerr spacetime. 



5 To convey an idea on the range of values that e could take, we use the results of Laarakkers & Poisson l'2(l for the quadrupole moment 
of rotating neutron stars. They show that the quadrupole moment could be approximated as q bj 2 . The factor b is a function of 
the mass and equation of state, with num eric al values varying from 2 to 12 (decreasing with mass and increasing with stiffness of the 
equation of state). It was emphasised in |20H that this relationship holds even for fast rotations. The fact that b is positive reflects 
oblateness of the mass distribution. Using this relationship we obtain: e = (6 — 1) j 2 , which tells us that for (stiff) rotating neutron star 
(with M = 1.0) e could be as large as 11 j 2 . 
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For reasons we discuss later in the paper, it is highly desirable for the quasi-Kerr metric to reduce to the familiar 
Kerr metric in Boyer-Lindquist (BL) coordinates in the limit e — > 0. The transformation between BL and the original 
H-T coordinates was provided by Hartle & Thorne 6 : 

i BL = t, 0bl = ((> (2.8) 

,-2i 



j_ 2 M 2 
~2r 

„'2 a ,r 2 



H3L = r- J - ir - r [(r + 2M)(r-M)-cos 2 e(r-2M)(r + 3M)} 



= r-j z M'f BL (r,8) (2.9) 
i 2 M 2 

#bl = 9- J — nr (r + 2M)sm6cose = - j 2 M 2 g BL (r,e) (2.10) 

Written in this manner, this transformation is not really practical for someone wishing to obtain the metric in BL 
coordinates. Fortunately, inverting the above equations with respect to the H-T coordinates is a trivial task. According 
to (|2.1U|) tbl = r + 0(j 2 ) and 9bl = 9 + 0(j 2 ); therefore we are free to use either set of coordinates in the functions 
Jbl t 9bl since these appear in the 0(j 2 ) terms. Thus, the inverse transformation is simply, 

r = r B L + j 2 M 2 fBh{rBh,0Bh) 

= #bl+J 2 A/ 2 .9bi>bl,0bl) (2.11) 

Hereafter, we shall drop the subscript on the BL coordinates. As the two sets of coordinates differ at order 0(j 2 ) 
the quadrupole-order pieces in 1)2. are immune with respect to the transformation l|2.11l) . Only the 'Schwarzschild' 
portion of the metric is affected by the coordinate transformation. For the contravariant metric components we obtain: 

Sht - .9^ + e (1 - 2A/A)- 1 [ / 3 (r) + / 4 (r) cos 2 9} 

+ e(l - 2M/r)[/ 3 (r) + / 4 (r) cos 2 6} 

[h 3 (r) + h 4 (r) cos 2 9] 

; r 2n [h 3 {r) + h 4 {r) cos 2 9} 
T l sm 

(2-12) 

where the functions /3, f±, /13, are given in Appendix 1X1 In eqns. (|2.12|) we have denoted as g^ 2 the 0(a 2 ) Kerr 
metric (in BL coordinates) with a = jM. From eqns. (|2.12|) one can extract the desired deviation h a p. 

Putting all the pieces together, the ansatz for building a quasi-Kerr metric in Boyer-Lindquist coordinates is: 



rr 
9ht 


— 5Ka 2 




- n" 


9wt 


- 5Ka 2 


4><t> 
Sht 


4>4> 


t<t> 
9ht 


— #Ka 2 



9ab = gfb + eh ab (2.13) 



where, 



h u = (l-2Af/r)- 1 [(l-3cos 2 6»)J : 'i(r)], h rr = (1 - 2M/r)[ (1 - 3cos 2 9) T^r)] (2.14) 

h 66 - -\[(l-3cos 2 9)T 2 (r)], h** = [ (1 - 3 cos 2 9) T 2 {r)] (2.15) 

r z sm 9 

= (2.16) 

The functions T\.2{r) are given in Appendix lAl 

Looking at the explicit form of h a p it is evident that these functions are divergent at r — * 2M. This would 
correspond to the location of the event horizon at the accuracy 0(j 2 ) of the H-T metric. Hence, there is no real 



Note that both the original transformation equations by Hartle & Thorne and their recent reproduction in Ref. |24| contain a sign error. 
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mystery behind this behaviour which is also present in the case of the full Kerr metric. However, since in our scheme 
we do not associate the quasi-Kerr metric with a black hole we are forced to be agnostic regarding the nature of the 
'surface' of the massive body. Whether is an event horizon or something else is an issue we simply do not address. 
This lack of information has no significant impact on the study of orbital motion of a test-body, as long as we do not 
consider orbits too close to 2M (assuming they exist). 

The situation is not so clear when one tries to study wave dynamics in a quasi-Kerr metric and compute (say) 
fluxes. To begin with, we should mention that a first approximation would be to completely neglect any gravitational 
flux impinging on the 'surface' assuming that is small compared to the fluxes at infinity. This is certainly true for 
the case of black holes [25|,[2(j and it is not unreasonable to expect the same in other scenarios. The main issue 
regards the computation of the fluxes at infinity, and their sensitivity to the absence of a concrete boundary condition 
at r ~ 2M. One could try different boundary conditions, say, perfect reflection or free propagation and assess how 
robust the flux results are (see for example Ref. ^(j). Other aspects of wave dynamics in a quasi-Kerr field, such as 
the existence and properties of quasi-normal modes, are extremely sensitive to boundary conditions and cannot be 
studied within our framework. 



III. GEODESIC MOTION IN THE QUASI-KERR METRIC 



A. The Hamilton-Jacobi equation: issues of separability 



An elegant method for deriving equations of motion for a point-particle in a given gravitational potential is the 
Hamilton-Jacobi (H-J) formalism (see Ref. [27j for background material and Appendix of the present paper). The 
end-product is the known H-J equation for the generating function S 1 , which in its general relativistic version is written 
as, 

1 ,8S_ dS_ 8S = 

2 9 dx-dx^ + dX [6 ' 

where A is the geodesic affine parameter. In essence, this equation follows from the definition of the Hamiltonian for 
point-particle motion, 

H(x a ,pp) = ^g^p ltPv = ~n a (3.2) 

by replacing the four momenta p a — dx a /d\ by p a — g a @ dS/dx 13 . Here the BL coordinates x a = {t,r,6,(f>} are 
conjugate to p a . 

One of the 'miracles' of the Kerr metric is that it allows full separability of the H-J equation, leading to the 
well-known Kerr geodesic equations of motion 0, ^5|. Separability with respect to X,t, (f> generates the conserved 
quantities /i (the test-body's mass), E (energy) and L z (angular momentum along the symmetry axis). This is always 
possible as long as the metric is stationary- axisymmetric. The special property of the Kerr metric (originating from 
its Petrov type- 1? character) is that it allows the additional, non-trivial, separability with respect to (r, 8) , leading to 
the third constant Q, the Carter constant. This is possible only in a restricted family of coordinate frames, the BL 
one among them [l7|. In fact, if we expand the Kerr metric with respect to a, then the H-J equation is separable at 
each individual order. 

On the other hand, the H-J equation is not fully separable with respect to the H-T metric (|2.6|l . not even when the 
Kerr limit e — > is taken. This reveals that the coordinate frame used in the original H-T metric is not 'privileged' 
and consequently a bad choice for a quasi-Kerr metric. Instead, BL coordinates seem like a good choice as in the 
limit e — > they admit separability. A similar conclusion can be reached for the general metric In its present 

form it does not allow separation of the H-J equation (or of the scalar and gravitational perturbation equations) even 
when the Kerr limit is taken (which corresponds to enforcing relation {EDO - 

The H-J equation for the full quasi-Kerr metric (|2.4|) is solved by assuming the standard form, 

S=^n 2 X-Et + L z <p + S r (r)+S e {9) (3.3) 

with an extra expansion, 

S r ,e = S r fi + eS r fi (3.4) 
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where S r ,e solves the exact Kerr H-J equation. After separating out the Kerr portion of (|3.1|) . and neglecting 0(ea, e 2 ) 
terms, we are left with, 



r (r - 2M) 



dr dr \ dr 



E 2 r*f 3 



where if is a separation constant and 

Z{r) = Kh 4 - U 



2M 



/i /•(/•- 2. V) ('f)' 



(3.5) 



(3.6) 



This non- vanishing function spoils the separability of (|3.5|l . In retrospect, this is not really surprising as the addition 
of the e-pieces 'contaminates' the Kerr metric in the sense that the resulting quasi-Kerr spacetime is not type-Z? 
anymore, although (in some loose sense) is arbitrarily 'close' (this is also the explanation for why the non type- 1? H-T 
metric does not lead to a separable H-J equation). 

At this point it would seem that no real advantage has being gained by using the quasi-Kerr metric (|2.4() as the H-J 
equation is still non-separable with respect to (r, ff). However, the fact that the equation is separable in the e — > 
limit can be exploited and by means of canonical perturbation theory |2?| ] can lead to the desired equations of motion 
(see Section lill B|) . 

Studying geodesic motion in a quasi-Kerr field is just the first step towards solving the EMRI problem. For 
gravitational wave observations it is pivotal to be able to compute accurate waveforms and fluxes. In the case of 
Kerr black holes, this is achieved by means of the Teukolsky equation ^(| (see Q, @ for reviews relevant to EMRI, 
and further references), which governs the dynamics of perturbations of the Weyl scalars {"00 > ipi} in the framework 
of the classic Newman-Penrose formalism. The remarkable result that all these perturbations are described by one 
single master equation is a consequence (once more) of the Petrov type-D character of the Kerr metric, which allows 
decoupling of the various perturbation equations 17]. Moreover, for a certain class of coordinate frames (the BL the 
most widely used one) the Teukolsky equation admits separation of variables, reducing the problem to the solution 
of simple ODEs. None of these properties survive for a non type-D metric such as (|l.lf) or Q2.4JI . However, as the 
quasi-Kerr metric is 'almost' type- D there may still be a chance of decoupling (at least partially) the perturbation 
equations for the Weyl scalars. This crucial issue is left for future work. In the unfortunate case that this expectation 
proves untrue, we may have to work with straight metric perturbations which are governed by ten coupled PDEs 
plus gauge conditions. These equations will still admit separation with respect to (f> so they can be studied in the 
time-domain in a '2 + 1' format. In the present paper, the issue of wave dynamics will be dealt by simply computing 
approximate waveforms 7 ( Section II V|) . but our conclusions should be valid even if we had rigorous waveforms at our 
disposal. 



B. Generic orbits 



Although we are mainly concerned with equatorial orbits in this paper, we nevertheless derive approximate equations 
of motion for generic quasi-Kerr orbits. Our objective is to demonstrate how one can arrive to such equations despite 
the non-separability of the H-J equation and without resorting to the full second-order equations of geodesic motion. 
This subsection is somewhat detached from the paper's main guideline and the reader should first consult Appendix FBI 
(where we review the action-angle formalism in the Kerr spacetime, following Ref. |28|) for definitions and notation 
used here. 

Our calculation consists of nothing more than the application of canonical perturbation theory (see |27| for further 
details) to point-particle motion in the Kerr spacetime. The full Hamiltonian can be written as the sum of the Kerr 
point-particle Hamiltonian amended with the quadrupolar perturbation, 

H = ~ gf Pa P0 + \e h a ? Pa P p = H + eH x , (3.7) 

The unperturbed Hamiltonian is a function of the actions only, H = H (Ip), while the perturbation depends on both 
actions and angles, Hi = Hi(w a , Ip). Expressing the Hamiltonian in this functional form is possible provided we 



A similar path was taken by Ryan who only calculated slow-motion, 'restricted Post-Newtonian', waveforms from circular equatorial 
orbits 
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have first obtained relations x K (w a ,Ip) and p K (w a , Ip). For the Kerr problem, these are given by implicit relations, 
see Appendix FBI 

In the unperturbed problem H the actions are constants and the angles vary linearly with respect to A. This is 
no longer true for the full perturbed Hamiltonian H. We seek a canonical transformation {w a ,Ip} => {ro",^} to a 
new set of angle-actions which would have exactly the above properties. Clearly, the old and new variables will differ 
at 0(e). The generating function takes the form [27j . 

F(w a J p ) = w k I k + ^(w a Jp) (3.8) 

Then we have, 

dF <9$ dF d$ 

I a = -P-=I a +-PL and ^ = ff_=^ + ±^. (3.9) 
dw a dw a dl a dl a 

Denoting as K(Ip, e) the new Hamiltonian we also have, 

^ = ™ => v" = P a \ + P<* (3.10) 
d\ dl a 

dL dK O + 0(e 2 ) (3.11) 



dX dw a 

where v a are the new fundamental frequencies. The function F solves the following H-J equation 

dF 

H{w a ,—)=K{I ,e) (3.12) 

The standard choice for the Hamiltonian K is, 

K(i p ,e) = H (i p )+e(H 1 )(i f3 ) (3.13) 
where we have defined the averaged perturbation with respect to the original angle variables, 

(ffi)(/0) = j>dw a H 1 {w k ,Ip) (3.14) 

Expanding the left-hand side of i|3.12[l we then find for the generating function, 

^^=em)-Ht] (3.15) 

The perturbed equations of motion are derived by first inverting eqns. (|3.9I) with respect to the old action-angles, 
isolating the secular changes in these relations and then substituting in the unperturbed equations w a (x k , Ip) (see 
eqns. (|B10|) in Appendix 151 ) . Note that in the function <!> we are free to interchange new with old actions at 0(e) 
accuracy. From l|3T5|l and the fact that &(w a ,Ip) is a periodic function of the angles w a [27| we find that, 

- 9® . , 

(3.16) 



'dw a 



which implies that there are no secular changes in the actions at 0(e). This also means that there is no secular change 
in {E,L Z , Q}. At leading order, quasi-Kerr orbits admit three integrals of motion, like in the Kerr geodesic motion. 

Averaging the perturbation in the angle variables (d<&/dl a (w k , Ip)) — (d$/dl a )(lp) ^ implies that in a secular 
sense w a = v A + (3 a where (3 a a constant phase. 

Therefore we have shown that the quasi-Kerr geodesic equations of motion can be directly derived from the cor- 
responding Kerr equations (|B10|I with the substitution v a v a in the left-hand sides only . Despite the identical 
shapes between Kerr and quasi-Kerr orbits, the difference in the frequencies will translate in a difference in precession 
rates for the periastron and for the orbital plane 8 . 



A basic feature of generic Kerr orbits is their characterisation by three incommensurate orbital frequencies Qi = Vi/vt (in BL time). 
This property results in orbits which are not closed, but rather exhibit precessional motions. Periastron precession is a consequence of 
O r ^ Q^, and orbital plane precession (Lense-Thirring precession) is a consequence of Qg ^ Clj,. Keplerian orbits, on the other hand, 
are fixed closed ellipses since their frequencies are degenerate, f2 r = Qg = fii. 
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Generic quasi-Kerr orbits will be discussed in detail in a future paper. Here we have shown that in a secular 
sense the motion has three conserved quantities (apart from fi), which is an important property and stands above 
of what one would expect for a general non-Kerr axisymmetric-stationary metric like ljl . As we mentioned in the 
Introduction, geodesic motion in such metric is not formally integrable and may well lead to chaotic behaviour in a 
test-body motion. 

In fact, orbital dynamics in general axisymmetric-stationary gravitational fields are typically much different than 
the ones in spherically-symmetric fields or even the Kerr field. We can easily demonstrate this difference, by asking 
the simple question: which axisymmetric-stationary gravitational fields can support circular orbits? In Newtonian 
theory the answer is that only a special family of potentials with the form V(r, 9) — V (r) + V±(0)/r 2 (where V , V\ 
arbitrary functions) admit orbits with u r — du r j dt = 0. In General Relativity the same requirement for a given metric 
g a b(r, 9) results in conditions involving the metric and its first derivatives. Once more, the Kerr metric is special as 
it satisfies these conditions. Further analysis of this issue will be presented in a future paper. 



C. Equatorial orbits 

For the special case of equatorial orbits the quasi-Kerr H-J equation JOJ is trivially separable. The resulting 
equation of motion are, 



r\2 



K) 



2 M , 
Vr(r) = (E 2 -1) + — - \L\ + a 2 (l - E 2 )]-^ + —(L Z - aE) 2 



1 + 2M , 



e l 



2M 



dX 

dt 
dX 



v 6 



F ' = A 



2M 



- (aE — L z 



E(r 2 + a 2 



2Ma 



(Ea - L z ) 



r J 



hE 



(3.17) 
(3.18) 
(3.19) 



In the standard manner, we express the test-body's radial location in terms of a pair of orbital elements, the semi-latus 
rectum p and eccentricity e, 



P 



1 + e cos x 



(3.20) 



In the weak-field limit (p 3> M) these parameters coincide with the familiar Keplerian parameters. The parameter 
X varies monotonically, while r 'runs' between the two radial turning points, the periastron r p = p/(l + e) and the 
apastron r a — p/(l — e). The elements (p, e) can be written as functions of E,L Z (and vice- versa) using V r (r a ) — 
V r (r p ) = 0. For the Kerr orbits V r (r) is a cubic polynomial and the third root r% (which is the smallest of the three) 
can be used to define a 'separatrix' of bound/unbound orbits. The transition occurs when r p = which defines 
an innermost stable bound orbit. In the quasi-Kerr case, the equation V r = is more complicated, but since the 
equations of motion (|3.19|) deviate only slightly from the corresponding Kerr equations, we expect the quasi-Kerr 
separatrix to be located close to the Kerr separatrix [2J| for the same {a,p, e}. 

Equatorial orbits are characterised by a pair of orbital frequencies Q r , f^. These are defined in terms of the radial 
period T r (the coordinate time required for the body to move from r p to r a and back), and of the periastron shift A(j) 
(total accumulated angle <f> over a period T r ), 



2tt 

J. , . 



A<j> 
TV 



(3.21) 



These frequencies are 'observables' as the emitted gravitational waveform (on an adiabatic timescale where backreac- 
tion is not significant) comes in discrete frequencies, 



(3.22) 



where k 7 m are integers (the latter associated with the azimuthal angle </>). The appearance of a spectrum of the form 
H3.22(l is a direct consequence of the periodicity of the motion. Note that in the Keplerian limit (p 3> M) we have 
A<p 27r, hence D, r — ► fK. 
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As a first stab at the problem of comparing Kerr vs quasi-Kerr EMRIs we calculate the above frequencies f2 r , 
for identical values {a,p, e} while varying the quadrupole deviation e. We use the simple expressions, 



T r = 



t n27T 

, U 

dt — = 



dx 



dr u 
dx u r 



(3.23) 



and 



A(j) 



dt 



27r , dr 

ax~, ~ 

dx u r 



2 71 



dx 



ep sin x 



vMx)) 



(i 



ecos 



(3.24) 



Results are shown in Figs. □ El & for the pair of orbits (p, e) = (10M, 0.5), (15 M, 0.5) and a = 0.5M. A quantity 
which nicely illustrates the difference between Kerr and quasi-Kerr orbits is the number of cycles Af required to 
accumulate tt/2 difference in periastron shift: 



Af - 



tt/2 



|A0 K - A0q K | 



(3.25) 



In Fig. |2| we show how Af varies with e and eccentricity e for fixed p and a. The dependence on e is quite weak, but 
there is a strong dependence on e. For example, for a moderate value |e| = 0.05 which corresponds to a fractional 
difference ~ 8% in the quadrupole moment, we only need about ~ 100 — 200 orbits to accumulate tt/2 difference in 
Acf>. Such large orbital dephasing would certainly manifest itself in the waveforms as we shall see in the following 
Section. 



i p=15M, e=0.5, a=0.5M 
i p=HM, e=0.5, a=0.5M 
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FIG. 1: Left panel: periastron shift difference A(/>k — A0 q K as a function of the deviation e. Right panel: number of cycles M 
required to accumulate 7r/2 difference in periastron shifts (defined by eqn. 13.251 ') as function of e. For both panels, we have 
considered two orbits: a — 0.5A/, e = 0.5 and p = 10M, p = 15M. 



IV. WAVEFORM COMPARISONS 



A. Kerr vs quasi-Kerr 

As we already pointed out in Section fill Al the calculation of rigorous gravitational waveforms and fluxes is a major 
challenge for any non-Kerr EMRI program. For the purposes of this study, which is just an initial investigation, it 
will suffice to use approximate waveforms, in particular the 'hybrid' waveforms discussed in |sl].|23l|. These waveforms 
are generated using flat spacetime wave formulae (see |l2j|.|29|| for more details) for the two components {h x ,h + }, 
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e 

FIG. 2: Difference in the radial period T r Kerr and quasi-Kerr orbits with a = 0.5M, e = 0.5 and p = 10M, p = 15M. 




FIG. 3: Number of cycles to accumulate tt/2 difference between Kerr and quasi-Kerr periastron shifts as a function of e and 
eccentricity e. We have fixed the other two parameters: p = 10M, a — 0.5M. 

coupled with exact relativistic geodesic motion. At least for the particular case of Kerr, it has been established that 
these waveforms agree quite well with rigorous Teukolsky-based waveforms |23|. Hence, it makes sense to assume that 
the same will be true for the quasi-Kerr case. 

In this Section we shall compare quasi-Kerr against Kerr waveforms with identical orbital parameters (p, e) , spin 
a/M and for the same kinematic initial conditions. Since a nonzero e imparts a change in the orbital frequencies we 
would expect this change to manifest itself as a significant cumulative phase difference between the corresponding 
waveforms after several orbits. This difference can be rigorously quantified in terms of the overlap function, that is, 
the scalar product between the normalised waveforms within the LISA sensitivity band: 



{huh2) = m rhmii df 

JO bhU) 



(4.1) 
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where a tilde denotes the Fourier transform of the waveform, a star stands for the complex conjugate, a hat denotes 
normalised waveform h(tf.) according to (hi, hi) — (h 2 , /12) = 1 and Sh(f) is the expected LISA sensitivity function 
(one-sided power spectral density). To calculate this latter function we use an analytic fit (discussed in Ref. (23|) to 
the numerically-generated sensitivity curve of Ref. |30l |. for one year of observation and omitting any confusion noise 9 . 
The overlap defined in l|4.1|l is closely related to the signal-to-noise ratio (SNR) [3l| , which is given by, 

SNR = (h, s) = A(h, h) = A = (s, s) 1/2 

Here we have assumed that the template h exactly matches the true signal s, and A is the amplitude of the signal 
in units of the normalised waveform h. An imperfect template or a disagreement between the template's parameters 
and those of the signal results in a reduced SNR, 



SNR = (hi, s) = A(h u h 2 ) < A 



(4.2) 



In this expression we have assumed that the normalised template hi does not exactly match the signal s = Ah 2 . In 
other words, a low overlap between a template and the expected signal implies a drop in SNR (which is accompanied 
by a drop in the event rate due to the decrease of the observable spatial volume) . 

Since we neglect any radiation reaction effect on the orbital motion our calculations will be consistent provided we 
'truncate' any waveform at the radiation reaction time scale Trr. We define Trr according to the following rule. For 
a test-body in Kerr spacetime we can include orbital backreaction using the approximate hybrid method § , j2^ , [12 . 
Then comparing Kerr waveforms with and without radiation reaction (denoted as /i^-(t), h RR (t), respectively) we can 
define Trr as the truncation time Tt r at which the overlap between these two waveforms drops below 95%. 



(h K (t 



2i 



RRJ 



hT (t = Trr) ) = 0.95 



(4.3) 



The particular threshold 0.95 was chosen as it coincides with the usual value of 'minimal match' for constructing 
template banks [U, [3j| . Based on our definition for Trr we can safely approximate the small body's motion as a 
geodesic for a time interval t < Trr. We have numerically computed Trr (see Fig. 01 as a function of fi/M, and not 
surprisingly, we find a linear dependence, Trr ~ M/[i. 
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FIG. 4: Comparing quasi-Kerr (solid curve) and Kerr (dashed curve) approximate hybrid waveforms for the orbit p — 10M, e = 
0.5, a = 0.5M and for e = 0.15. In addition, we include a Kerr waveform taking into account backreaction on the orbit 
(assuming the same initial orbit), represented by the dashed curve. All waveforms are shown at a time window close to the 
radiation reaction timescale Trr. 



We have observed that addition of confusion noise does not affect the overlaps, although it slightly affects the estimated radiation 
reaction time scale. 
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For the purpose of waveform comparison we select the moderately relativistic orbit p = 10M, e = 0.5, a = 0.5. For 
a mass ratio n/M ~ 10~ 5 we have T KK w 27500A/. The quadrupolar deviation is fixed at e = 0.15 which corresponds 
to a fractional difference SQ/Q K = 40%. The actual waveforms are shown in Fig. 0] Clearly, after a time lapse 
~ Irr the accumulated phase-difference between the two waveforms is quite significant. The waveform deviation due 
to a nonzero e can be also compared against the deviation imparted by allowing the orbit to evolve under radiation 
reaction. For the particular case of Fig. after a time interval t ~ Trr, the waveforms dephase mainly due to the 
non-Kerr property, rather than due to radiation reaction. 

In order to quantify the waveform difference in a more rigorous manner we compute overlaps between Kerr and 
quasi-Kerr waveforms truncated at t = Trr. Some results are given in Fig. [5] as a function of the mass ratio for the 
two values e = 0.07 and e = 0.15. One can see that the overlap can drop quite dramatically (down to 40%) even for 
small values of e (0.07) and for realistic mass ratios (~ 10~ 6 ). 

This an important result of the present paper: The typically low overlaps between quasi-Kerr and Kerr waveforms 
(for the same orbital parameters and spin) simply mean that using Kerr templates for studying EMRIs in a true 
non-Kerr spacetime might result in low SNRs, accompanied by significant loss in the number of observed events. 




log 1() (u/M) log ]0 (u/M) 

FIG. 5: Right panel: radiation reaction timescale Trr (see main text for definition), as a function of mass ratio pt/M. Left 
panel: overlaps (expressed in %) between quasi-Kerr and Kerr waveforms for the same orbit p = 10M, e = 0.5, a = 0.5M and 
for e = 0.07, 0.15, truncated at t — Trr(/i/M). The solid line is the spline interpolation between the data points. 



B. Prelude to the 'confusion problem' 



According to our results, it is clear that even a modest (say ~ 10%) deviation from the Kerr quadrupole moment is 
sufficient to render any Kerr template waveform inaccurate assuming the same orbit {p, e} and spin parameter a/M. 

This statement begs the question: is it actually best to compare waveforms corresponding to the same set of 
parameters {a,p, e} ? One could argue that we should have rather fixed the orbital frequencies, 



nf = nf, i = ^,r (4.4) 

Then one would expect a good agreement in phase between the two waveforms. Note that this choice would correspond 
in comparing waveforms with different {p, e}, given fixed values of a, e. If (p, e) and (p, e) are the Kerr and quasi-Kerr 
orbital elements, respectively, and assuming a small deviation Sp — p — p ~ O(e) (and similarly for e) we have, 

nf ( P , e ) = nf (p, g) = nf (p, g) + e siW ( P , e ) 

^ip+^fe=- £ flW(p ) e) (4.5) 
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This 2x2 system has a non- vanishing determinant and a unique solution for (Sp, Se). Note that for e = this solution 
becomes trivial, but this is true as long as we do not allow for variation in a/M. In Fig. [5] we show the solution of 
(j4.5(l for a = OAM, e — 0.1. We can see that the most pronounced deviation between Kerr and quasi-Kerr orbits 
occurs in the small p, small e region. This is understandable as in this case the body spends most of its orbital time 
in the strong-field where the deviation from the Kerr metric is stronger. As expected for p 3> M both Sp, Se — > 0. 

Our preliminary results on waveform comparison are quite alarming: fixing the orbital frequencies instead of {p, e} 
lead to overlaps very close to unity between Kerr and quasi-Kerr waveforms. For the particular case of n/M = 10 -5 , 
a = 0.3M and e = 0.1 and for {p, £} = {TOM, 0.3} we get from gSJ {p,e} K = {9.906M, 0.317} and the resulting 
overlap is 97.6% (for the time interval Trr w 31000M). Hence, we have naturally stumbled upon a case of the 
confusion problem: the possibility that for a given non-Kerr EMRI waveform there might be a corresponding Kerr 
waveform with a different set of parameters which could mimic the former (that is produce an overlap > 95%) on a 
radiation reaction timescale. 

Here we have given one example where this scenario is indeed true, but this is not the end of the story. Allowing a 
time window, longer than Trr, for which radiation reaction effects become important the overlap could significantly 
reduce as a consequence of different orbital evolution between Kerr and quasi-Kerr. Orbits which are initially 'tuned' 
to have almost identical orbital frequencies are not expected to preserve this property for t > Trr. 

The non-trivial issue of waveform confusion is of great importance for LISA and certainly requires further investi- 
gation. We will address it in more detail in a following paper. 
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FIG. 6: Deviation in the orbital elements p, e between Kerr (squares) and quasi-Kerr (circles), after equating the orbital 
frequencies Q r , Q^, and for fixed a = OAM, e = 0.1. 



V. CONCLUDING DISCUSSION 



One of the LISA mission main science goals will be the 'mapping' of Kerr spacetime. In order to perform this 
crucial test of General Relativity we should be in a position to prepare data analysis tools capable of gauging any 
possible deviations from the Kerr metric. In this paper we have introduced a practical scheme which could be a strong 
candidate for this task. Our scheme is based on the notion of a 'quasi-Kerr' metric, that is, a metric that deviates 
only slightly from the known Kerr metric, while still being stationary and axisymmetric. This deviation can first 
appear in the value of the spacetime's quadrupole moment. We only consider such deviations and ignore any other 
possible disagreement in higher multipole moments. 

Our quasi-Kerr metric is built with the help of the exterior Hartle-Thorne metric, which describes the spacetime 
outside any slowly-rotating, stationary axisymmetric body. Given this metric, we have first studied geodesic motion 
of a test-body. Despite the fact that the Hamilton- Jacobi equation is non-separable, we manage to derive approximate 
equations for generic orbits employing canonical perturbation theory. 
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In this paper we mostly focused on equatorial orbits (in which case the Hamilton-Jacobi equation is trivially 
separable) and for given orbital parameters and spin we computed the change imparted to the orbital frequencies 
and periastron advance by the non-Kerr metric deviation. We found that for (say) a ~ 10% fractional difference in 
the quadrupole moment, the cumulative change in these quantities is significant after about ~ 100 orbits, and even 
for moderately relativistic orbits, it significantly affects the gravitational wave phasing. This was demonstrated by 
generating approximate 'hybrid' waveforms (for orbits with identical p, e), for both Kerr and quasi-Kerr spacetimes, 
and computing their mutual overlap function. 

These overlaps typically drop below ~ 60% for the relevant mass ratio (allowing for a time interval at which any 
radiation reaction effects are still not important). Such low overlaps mean that we might expect significant loss in 
signal-to-noise ratio if we attempt to use Kerr waveform templates for the detection of a gravitational wave signal 
generated by an EMRI in a quasi-Kerr field. This statement holds as long as we compare waveforms with the same 
orbital elements. However, the picture becomes more blurred when we compare waveforms which correspond to 
different orbits but same orbital frequencies. Then it is possible for a pair of Kerr and quasi-Kerr waveforms to match 
extremely well. This confusion problem could be a generic feature of any non-Kerr data analysis effort, and it means 
that we may confuse a possible true non-Kerr event with a Kerr EMRI with the wrong parameters. A possible way to 
break this degeneracy is to include backreaction effects, in which case orbits with initial identical frequencies diverge 
on the radiation-reaction timescale. This issue will be further analysed in a following paper. Similarly, we shall extend 
the present calculations to the case of non-equatorial orbits, but we expect that the conclusions of the present paper 
will remain true in this more general scenario. 

The scheme we propose here clearly requires further development, especially towards the direction of computing 
rigorous quasi-Kerr inspiral waveforms. The subject of wave-dynamics in a general vacuum axisymmetric-stationary 
spacetime or even a less general spacetime like the quasi-Kerr, is almost unexplored. Work in progress aims to inves- 
tigate if it would be useful to work within the Newman-Penrose formalism |l5j , trying to formulate some approximate 
'quasi- Teukolsky' equation, or if it would be more practical to work directly with metric perturbations. 
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APPENDIX A 

Here we provide the explicit form of the various functions appearing in the H-T metric, eqn. I|2.6[) . 
Fi{r,0) = - 8Mr t(l_ 2M) [ (r - M)(16M 5 + 8M 4 r - lOAfV - 30Afr 4 + 15r 5 ) 

+ cos 2 8(48M 6 - 8M 5 r - 24AfV - 30Af 3 r 3 - 60A/V + 135A/r 5 - 45r 6 )] + A x (r, 9) 

= /i(r) + / 2 (r) cos 2 0, (Al) 

F 2 (r,0) = _5(,-Af)(2Af 2 + 6M r 3, 2 ) Qg2 
2V ' ; 8Mr(r-2M) X 

= h(r) + h(r) cos 2 9, (A2) 
n i m (P - 72Af 5 r) - 3 cos 2 6»(L - 56A/ 5 r) . . n . .. . . „ „ 

Gi(r,0) = ^ 8Afr 4 (r-2Af) = 5i W + fla(r) cos 2 , (A3) 

„, as (16Af 5 + 8Af 4 r-10A/ 2 r 3 + 15A/r 4 + 15r 5 ) M „ 2/11 , . ... 
#iM) = >-{l - 3 cos 2 9} + A 2 (r, 9) 

= hi(r) + h 2 (r) cos 2 9, (A4) 
5f2M 2 — 3Mr — 3r 2 ) 

H 2 (r,9) = J — ^{l-3cos 2 0}-A 2 M) = ft 3 (r) + fc 4 (r) cos 2 9, (A5) 



where 



P{r) = 15 r b - 45 Mr 5 + 20 A/ 2 r 4 + 10 Af V + 8 Af r 2 + 80 AP (A6) 
Mr,0) = 15r i~l M) ln(^—]{l-3cos 2 8}, (A7) 



16M 2 \r- 2M 

Mr,0) = ~ 15(r l~ 1" ( —777 ) ! ' - ■> »^ ^ • (AM 



16Af 2 V r -2M 
The radial functions -Fi^r) appearing in the quasi-Kerr metric (|2.13|) are: 

*(r) = - 5(r r M L (2M 2 + 6Mr - 3. 2 ) - i^Ml ln f _!_) (A9 ) 
W 8Air(r - 2Af ) V y 16M 2 ^ - 2M/ 

jr ( r ) = _JL (2Af2 _ 3 Mr - 3r 2 ) + — ^(r 2 - 2M 2 ) ln I ) (A10) 

y ' 8Mr y ' 16Af 2V ' \r-2M) y ' 

APPENDIX B 

In this Appendix we p resent, in some detail, the Hamilton-Jacoby/action-angle theory for a test-body in the Kerr 
spacetime, see Refs. [IJ, [2^, [34| for full exposition to the topic. 

The starting point is the Hamiltonian for point-particle motion (a, j3, k € {t, r, 9, </>}), 

Ho(x a ,pp) = -g^PaPp = -T^ 2 (Bl) 

Generating a canonical transformation {x a ,pp\ => {Q a ,jfi}, with the requirement of a vanishing new Hamiltonian 
K, leads to the following H-J equation for the generating function S(x a ,jp, A): 

1 dS dS OS 

2 9 d^d^ + d\~° (B2) 
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As we have emphasized, the Kerr metric allows full separability of this equation in BL coordinates. We write 

2 2 

S(x a , 7/3, A) = !jX - Et + L z + S r (r) + S e (9) = fjX + W(x a , 7/3 ) (B3) 
The function W(x a ) <jp) itself is a generator of a canonical transformation with K = — fi 2 /2. In any case we find, 

S r (r) = J dr ^ and S e (9) = J deVo (B4) 



where A = r 2 — 2Mr + a 2 . 
well-known expressions (see 



The potentials R(r,ak) and 0(0, afe), where afe = {— /i 2 /2, i?, L z , Q} are given by the 



ine pc 

MM 



R= [(r 2 +a 2 )E-aL z ] - A [fi 2 r 2 + (L z — a E) 2 + Q] 



9 = Q-cos 2 



a 2 (fj 2 -E 2 ) + - 



(B5) 



The 'third' constant Q is the familiar Carter constant and is related to the original (r, 9) separation constant C as 
Q=C 2 -L 2 z -a 2 E 2 . 

Equations of geodesic motion are derived from, 

pt = -g tt E + g t ' t >L z = (r 2 +a 2 )A- 1 T-a{aEsm 2 9-L z ) 

p r = ± g rr A -l y/R = ±Y,~^ y/R 

p 4> = g 't"t>L z -Eg t ^ = a A" 1 T-aE + L 2 z / sin 2 9 



a n a $ 



dS 



where T = (r 2 



i E — a L z and S = 



(B6) 

a 2 cos 2 9. Similarly, the covariant momentum components are: pt = 



—E, p r — ±\/i?/A, po — ±v6, Pcj, — L z . Incidentally, we point out that an equivalent 'integrated' form of the 
equations of motion can be written by using Q a = dS/d"f a , where we can always make the identification 7^ — a^. 

The periodic character of the motion is revealed when we introduce action-angle canonical variables. The actions 
are defined as (hereafter i £ {r, 9, </>}), 



Ji = <p pidq 1 => < 



f J r — 2 dr sfR/A 
J e = 2 J®' d9 Ve 

k Jcf, — 2 7T L z 

and are constant functions Ji(a K ). In the expression for Jg we have used 9 ni 9 S for the turning points of the 9- motion. 

The next step is to generate the canonical transformation {x a ^pp} =>■ {w a , Ip} with Ip = {p t , J r , Je, J^}- Clearly, 
Ipiotn), or a K {Ip). Since the new momenta J3 are constants with respect to A the new Hamiltonian must be K — 
K{Ip) = K(ak). Hence we can set K — — /i 2 /2 and use W(x a ,I/3) of (|B3|) as the generating function. 

For the new coordinates conjugate to Ip we have, 



dw a 
~d\ 



ok _ aHo 

dl a dl a 



(B7) 



with v a ,P a constants. As expected, w a are angle variables and v a fundamental frequencies. These frequencies are 
related to the orbital frequencies (in BL time) as, 



Mfli = Vi/vt 

The explicit expressions can be found in Ref. |2Sj. We also have, 

a _dW _dW da K 
dl a da K dl a 



(B8) 



(B9) 
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Expanding eqns. I|B9 



w 



u t X + f3 t 



-2u t 



w r = v r A + [3 r = —2 v r 



ve A + p 6 



-2 v$ 



8S r 


dS e ' 


dpt 2 


dp 2 


~dS r 


dS e ~ 




dp 1 


~dS r 


dSg 


dp? 


dp 2 


~dS r 


+ dSe_ 


dp, 2 


dp 2 



+ t 



dQ 
dJ r 

dQ 
dJ e 

dQ 

dJi, 



dS r dSg 
~dE + ~dE 

dS r dSg ' 
~dQ + ~dQ_ 

dS r dSg ' 
~dQ + ~dQ_ 

dS r dSg 

~d~Q + ~d~Q 



1 
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dS r 

dL z 



dSg 
dL z 



(BIO) 



These are the 'integrated' form of the equations of motion and they also provide the x a (w k , Ip) relations in an implicit 
form. Eqns. IjBlOp are fully equivalent to eqns. JB6J); this can be easily demonstrated by differentiating (say) the w r 
and w e equations, 



r- iT 1/2 r + a 2 cos 2 9" 1/2 6 + - n r 



r 2 R-^r + ^^ee-^e + ^ng 



Q- 1 ' 2 o - r- 1 ' 2 e 



where here an overdot stands for d/dX and 



1_ dQ_ 

Vi dJi 



Combining eqns. (|B11|) . 

(Kr - Kg) 



8- 1 / 2 - R- 1 ' 2 f\ = Q- 1 ' 2 9 = Rr 1 ' 2 f = a(r, . 



(Bll) 



(B12) 



(B13) 



since n r ^ Kg for non-equatorial orbits. The remaining function a{r, 9) is subsequently specified by substitution into 
either of eqns. (|B11I) . We find a = (r 2 + a 2 cos 2 ff)^ 1 = Hence, we have arrived at the expected (r, 9) equations 
of motion. The remaining (t, <j>) equations can be obtained in a similar fashion from the (w t , w^) equations. 
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